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Abstract 

We investigate a two-parameter entropy introduced by Schwammle and Tsallis and obtain its 
probability distribution in the canonical ensemble. The probability distribution is given in terms 
of the Lambert W-function which has been used in many branches of physics, especially in fractal 
structures. Also, extensivity of Sq^qi is discussed and a relationship is found to exist between the 
probabilities of a composite system and its subsystems so that the two-parameter entropy, Sq^qi, is 
extensive. 

1 Introduction 

It is known that some physical systems cannot be described by the Boltzmann-Gibbs(BG) statistical 
mechanics. Within a long list showing power-law behaviours, we may mention diffusion [1], 
turbulence [2], transverse momentum distribution of hadron jets in e"'"e~ collisions [3], thermalization 
of heavy quarks in collisional process [1], astrophysics [5], solar neutrinos [6], and among others 
[71 [HI [9l [ini [H]. Such systems typically have long-range interactions, long-time memory, and 
multifractal or hierarchical structures. To overcome at least some of these difficulties, Tsallis 
proposed a generalized entropic form \\.2\ [T3l [T^ , namely 

Sq = k ^^=^^' , (1) 

where, /c is a positive constant and uo is the total number of microscopic states. It is clear that 
the q-entropy {Sq) recovers the usual BG-entropy {Sbg = —kJ2i'=iPi^^Pi) limit q ^ 1. By 

defining q-logarithm 

lug X = — (hii X = In x) , (2) 
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the entropic form, Sq, can be written as 

Hence, in the case of equiprobability, Pi = the well-known Boltzmann law is recovered in the limit 
q —f 1. The q-entropy also satisfies the relevant properties of entropy like expansibility, composability, 
Lesche-stability [T7] and concavity (for q > 0). The inverse function of the q- logarithm is called q- 
exponential [1^ and is given by 

1 

expg X = [1 + (1 — q)x] 1-9 (exp]^ x = exp x) . (4) 

Recently in [18], the two-parameter logarithm liiq^qr{x) and exponential exp^ were defined which 
recovered q- logarithm and q-exponential, in the limit (; — > 1 or — > 1. So, the two-parameter 
entropy, similar to Eq. ([T]), can be defined as 



Sq^q' 



^pan,,,,- = ^— -^Pi[exp(y-^(p^-^-l)) -l] . (5) 



i=i Pi - -i 



The above entropy for the whole range of the space parameter does not fulfill all the necessary 
properties of a physical entropy. Therefore, it may not be appropriate for describing the physical 
systems, but useful for solving optimization problems. In this paper, we will find the probability 
distribution pi for the two-parameter entropy Sq^qi (jl7p . when canonical constraints are imposed on 
the system. As a result, it will be shown that the probability distribution is expressed in terms of 
the Lambert W- function |19 t I2UI I22j. also called omega function, which is an analytical function 
of z defined over the hole complex z-plane, as the inverse function of 

z = We^ . (6) 

Using this function, it is possible to write the series of infinite exponents in a closed form [23] 

■ = -"-•<-'■'-' . (7) 

mz 

The above equation shows that Ty-function can be used to express self-similarity in some fractal 
structures. This has been recently shown in a number of studies. Banwell and Jayakumar [23] 
showed that Vl^-function describes the relation between voltage, current and resistance in a diode, 
Packel and Yuen [25] applied the TF-function to a ballistic projectile in the presence of air resistance. 
Other applications of the W- function include those in statistical mechanics, quantum chemistry, 
combinatorics, enzyme kinetics, vision physiology, engineering of thin films, hydrology, and the 
analysis of algorithms [271 128j . Sq may become extensive in cases where there are correlations 
between subsystems. Hence, a point to be discussed will be the extensivity of Sq^qi [25]. 
This paper is organized as follows. In Sec. 2, we will review the Generalized two parameter entropy, 
Sq^qi and its properties. In Sec. 3, the probability distribution of the two-parameter entropy, Sq^qi, 
will be obtained in the canonical formalism. In Sec. 4, assuming that the two-parameter entropy be 
extensive, we will develop a relationship holding between probability in the composite system and 
probabilities of subsystems and in Sec. 5, we will have a conclusion. 
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2 Generalized two parameter entropy, 5^ g/ 

In this section, we will review the procedure for finding the two-parameter entropy, Sq^qf [18j . As 
we know, it is possible to define two composition laws, the generalized q-sum and q-product |16j . 
defined as follows 

X ®qy = X + y + {I - q)xy {x ®i y = x + y) , (8) 

X (^q y = {x^^'^ + y^^'^ - 1)^ {x®iy = xy). (9) 

Using the definition of q-logarithm (Eq. ^) and q-exponential (Eq. ([3])), the above relations can be 
rewritten as 

ln(;(a;y) = In^x In^y , (10) 

\nq{x (^qy) = \\lqX + hlqy . (11) 

Very recently in [18j, Eqs. (jlOp and (jlip were generalized by defining a two-parameter logarithmic 
function, denoted by lug^g' x, which satisfies the equation 

lng,g/ (X ®q y) = lUg^g/ X ®ql lUg^g/ ^ . (l2) 

Assuming In^^g' x = ^(Ingj;) = g[z) and using x = y u\ Eq. (fT2l) . then applying some general 
properties of a logarithm function like 

lng,5' 1 = 0, (13) 

■^lUq^q: X\x=l = I , (14) 

the two-parameter generalized logarithmic function will be given as 

exp ' ' 



(^(xi-" - 1)) - l] = In,, e'^"^ , (15) 



and its inverse function will be defined as a two-parameter generalized exponential, 



1 



exp,,,,x = {l + ^— 4ln[l + (l-g')^]}"' . (16) 



The entropy can be constructed based on the two-parameter generaUzation of the standard logarithm 

^ \ k ^ 1 ' 

Sq,q' =kY.Pi\nq^q' - = 3 II '^''P ( 1 '^^^^^^^ ~^^) ~ ' ^^''^ 

j=l Pi ^ ^ i=l ^ ^ Q 

which, in the case of equiprobability {pi = ^ Sq^q' = kluq^qiuj. 

The above entropy is Lesche-stable and some properties such as expansibility and concavity are 
satisfied if certain restrictions are imposed on {q,q'). 
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3 Finding probability distribution in the canonical ensemble 

In this section, we are interested in maximizing the entropy Sq^qi under the constraints 

Y.m-i = Q, (18) 

i=l 

iu 

Y,p,e^-E = {). (19) 



i=l 



These constraints are added to the entropy with Lagrange multiphers to construct the entropic 
functional 



Sq,q' 



k 



+ Q(^K-l)+/3^Pi(ei-S) = 



(20) 



i=l i=l 

To reach the equihbrium state, the entropic functional $<j,g' should be maximized, namely 



dpi 



=^ lng,g/ — +pi ^ 

Pi op- 



1 «'""'(i)+„+^(,._E)=„, 



(21) 



Using the definition of the two-parameter logarithm (Eq. (jlSp ) and after some calculations, we get 
{^Y^W^ - 1)) (l - (1 - Q')P^''^^) = 1 - (1 - g') (« + die, - E)) . (22) 



exp 

To solve the above equation and to find pj(ej), one may define 

1 



l-q 

and so, Eq. ([22]) can be rewritten as 



(K^-^-l) ^ l-(l-g>,'?-i = (g-l)z, + g', 



{q' +{q- l)zi) exp(zi) = ji , 



with the definition 



ji = l-{l-q')(^a + (3{e^-E)) . 
Solving the above equation gives us 



Zi = W 



re 9-174 



q-1 



+ 



1 



(23) 



(24) 



(25) 



(26) 



where, W{z) is the Lambert VF-function [H [201 [2T1 122] . From Eqs. ([26]) and ([23]), the probability 
distribution is given by: 



R = i-{l + (l-,)M'[^^]}"' 

Z „ I- L cr - 1 J J 



q-1 



where, 



g 9-1 

[-6-3-1 7j 



^. = E{i + (i-^)^[^]} 
i=i ^ 



(27) 



(28) 
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Figure 1: Normalized probability distributions are depicted for different values of q and q in (a) and (b), respectively. 
It is also shown that for q = 1.1 and q = 0.9, the plots approach to normalized Guassian. 

In Eq. psp . /3 is entered as an inverse of pseudo-temperature, but it may be interesting to write the 
probability distribution in terms of a deformed q-exponential which is more similar to the Boltzmann 
probability distribution. So, 7 in Eq. (22) can be written as 



7. = (1 - a(l - .0) (1 - T^^(^(-^ - E)) - (1 - a(l - 
where, f3q may be defined as the inverse of the pseudo-temperature 

/3g 



exp (-/3g(ei - E)) 



1 



ksTq 1 - a(l - q') 
The probability distribution can be written in a better form 



(29) 



(30) 



w 



69-1(1 - a(l - q')) 



1 



(expg(-/35(ei - E)) 



l-q 



(31) 



We can assume the energy level, quadratic function of the variable Xj. The continuous 

normalized probability distribution as a function of x can then be rewritten as: 



p{x) 



{1 + (1 - q)w[^{l - (1 - gO(a + P{x' - E)))]}^-' 



IZo {l + {l-q)W (1 - (1 - q'){a + P{x^ - E)))]} dx 



(32) 



This is illustrated in Fig. 1. Comparing these plots with the normalized Guassian distribution shows 
that in the limit q,q' 1, the plots approximate the Guassian function as expected. 
It is also possible to repeat the procedure with the energy constraint J2^iPi^i — E. But in that 
case, the equation which maximizes the entropy, is not solvable. 



4 How to interpret the entropy Sq^q' extensive? 

Extensivity, together with concavity, Lesche stability, and finiteness of the entropy production per 
time, increases the suitability of an entropy. But, is BG-entropy the only extensive one? In Ref. |29j . 
a question is raised as to whether entropy Sq is extensive or not? The answer is meaningful only if 
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(a) (b) 

Table 1: (a) shows the independence of two equal two-state subsystems A and B. The joint probabihties are given 
by the multiphcation of subsystems, probabihties. (b) shows an example of the correlated A and B. 

a composition law is specified, otherwise, it is tacitly assumed that subsystems are independent. 
Special correlations can be mathematically constructed such that Sq becomes extensive for an 
adequate value of g ^ 1. For example, we can consider two equal subsystems A and B with 
probabilities pf and , respectively. The probabilities of a composite system can be described by 
joint probabilities pfj^^ ■ In Tables (a) and (b), two cases are considered: (a) shows the independence 
of subsystems A and B, namely Pif~^ = pfpf and thus, the BG-entropy and the q-entropy become 
extensive and non-extensive, respectively. 

Sbg{A + B) = Sbg{A) + Sbg{B) , (33) 
Sq{A + B) = Sq{A) + Sq{B) + (1 - q)Sq{A)Sq{B) , (34) 

However, (b) shows a special correlation between subsystems which leads to non-extensivity for the 
BG-entropy and extensivity for the Tsallis entropy at g = 0. 

So{A + B) = So{A) + So{B) . (35) 

As can be seen in Table, (b), one of the states of the composite system appears with the 
zero probability and so, the number of effective states is w'^f'f^ = 3, which is not equal to 
w^~^^ = X = 4. This simple model can be improved to describe non-ergodic systems, 
where not all the states are accessible. In the following, along the lines of what is done in [291 130|. 
we will find a relation between the probabilities of the composite system (joint probabilities) and 
the probabilities of subsystems (marginal probabilities), such that the two-parameter entropy Sq^qi 
is extensive. 

Now, we are interested in finding the extensivity condition for the two-parameter entropy. Consider 
N subsystems {Ai, A2., ■ ■ ■ ,^Ar), each with the probability pi^ (s is related to that system). The 
probabilities in the composite system, pf^ti^'^t^'^'^'^ ■> should satisfy the condition 

The marginal probability related to the system, s, is defined as 

Pis - Ph,i2,-,iN ■ y^'l 

ii,i2,...is-iis+i---,iN 

It may be interesting to find the condition which makes the entropy Sq^qi extensive. In other words, 
we want to know the relationship between the probability in the composite system and probabilities 
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of the subsystems when the entropy Sg^q' is extensive. Let us consider the relation 

1 ^1 

A,+A2+...+A^ = ^^P?,?' ( E 1%.?' -47 + '/'n,i2,...,*iv) > (38) 

where, 4>ii,i2,...,ipf is set to ensure Eq. ()36p . The above equation can be rewritten in a different form 

1 11 1 / ^ / N 

Ai+A2+...+A^ = — ^9.'/' — • • • -Aj; ®q,q' exp,,q. ((/>n ,^2,...,^^^ ) , (39) 

Pil,i2,...,iN Pil PiN 



with the definition of (8'g,g'-product [18] 

X ®q,q' y = eWq,q' (lng,q' X + lUg^g' y) . (40) 

A nonzero function (pi-i^, 12,..., ipf is related to the existence of the correlation in the system, because in 
the case of independent subsystems {q,q' —i- 1), the (8)g,g'-product becomes the usual product and 
021,42, ■••,ijv ~ 0. In Ref. [30], the values of (pij, making the entropy 5^ extensive, are obtained for two 
equal two-state subsystems. In the case of equiprobability, Eq. ([39|) save for the function (pii^i2,...,ij^, 
shows the generalized multiplication of the number of states of the subsystems, which may be defined 
as the effective number of states 

W^^f = W^^ W^^ (g)g,,/ . . . (E)q,q' W^^ Og,q/ eXPg^g, ((/^j, ) , (41) 

The entropy of a composite system similar to Eq. (jl7p can be defined as follows 

N . 

Sq,q{j:^s)^k E Pt::^:^-^^''^..'^,TAkTA^- (42) 

S=l li,J2,...,«jv P«l,l2,...,«JV 

Using Eq. (p8|) . the entropy can be written as 

AT N 



S=l 11,12, ...jiN s=l Pis 

N . 

- 2^ Ph,i2,-,iN 2^^%?' As ^ 2^ Pii,i 



N 

AI+A2+...+AN 

_ _ ^,i2,--;iN Yll,l2, — ,iN 

1 Pi. 



ll,l2,...,tN S=l ^ts Ji,l2,...,«jv 

N -1 

^EE^^^K,'47 + ^ E Pttt't.■^^''<t>^.,n,...,^. 

s=l is Pis ii,i2,...,iN 

N 

Y^Sq,AAs) + k l'tw;-'"^"<^n,.2,...,.. , (43) 

s=l ii,i2,...,ijv 



where, the definition of marginal probability Eq. (j37j) is used in the last line. Eq. (j43]) ensures 
extensivity of Sq^q' if the constraint 

^Ai+A2+...+Ajv, 

ll,i2,...,iJV 

is satisfied. In other words, assuming the above constraint will be equivalent to the existence of 
extensivity. According to Eqs. (jlSp . ()16p . and (I38p . for the probability of a composite system, we 



E Pn:.w;-^^"'^n,.2,...,2. = , (44) 
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get 



A1+A2+...+AN 

Pil,i2,...,iN 



1 

{1 + ^ In [1 - AT + (1 _ - E ^Mj^ipV - 1)]] } • (45) 



It is clear that in the hmit (? — > 1, the proposed relation of probabilities is recovered [291 [30] 



A1+A2+...+AN 

Pii,i2,---,iN 



1 

TV q'-l 



(46) 



where, in the limit — > 1, the usual product of probabilities, Pi-^i^' 
which describes the case of independent subsystems. 



s=l 



n^iPi., is given, 



5 Conclusion 

In this paper, a special set of two-parameter entropies [TH] were maximized in the canonical ensemble 
by the energy constraint Yl^=i Pi^i = We expected that the probability distribution, Pi{ei), can be 
expressed in terms of the generalized two-parameter exponential defined in [18]. But unexpectedly, 
solution of the related equation took the form of the Lambert function which has been used in many 
branches including statistical mechanics, quantum chemistry, enzyme kinetics, and thin films, among 
others. The Lambert function can also describe some fractal structures because infinite exponents 
can be written in terms of the Lambert function (Eq. [7]). This suggests that the two-parameter 
entropies are probably related to the fractal structures in a phase space which may be a subject 
for future study. Also, assuming extensive Sg^qi, the probability of a composite system was given in 
terms of probabilities of the subsystems. 
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